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The  stress  distribution  in  bond  layers 
of  two  different  thicknesses  (50  pm  and 
200  pm)  was  calculated  by  finite  ele- 
ment analysis  for  pairs  of  rectangular 
cross  section  metal  bars  bonded  to  each 
other  and  subjected  to  four  point  bend- 
ing. These  stresses  were  used  to  aid  in 
identification  of  the  failure  origin  by 
use  of  the  Weibull  risk-of-rupture  (RR) 
function.  By  use  of  the  stress  distribu- 
tions, the  characteristic  strength  from 
50  pm  bond  test  specimens  could  be 
correlated  with  that  for  200  pm  bond 
test  specimens  when  the  failure  was  as- 
sumed to  have  an  interfacial  origin. 

The  finite  element  meshes  were  refined 
twice  and  the  ratios  of  characteristic 
strengths  were  recalculated  and  re- 
mained virtually  unchanged  by  each  of 
the  mesh  refinements.  Hence,  the  iden- 
tification of  the  interface  as  the  failure 


origin  remained  consistent.  Further,  the 
use  of  stresses  extrapolated  to  zero 
mesh  size  also  produced  the  same  ra- 
tios. Therefore,  the  RR  calculations  do 
not  appear  to  be  sensitive  to  the  mesh 
sizes  used  for  the  stress  calculations 
when  the  meshes  are  comparable  or 
when  changed  in  a comparable  manner. 
The  results  show  this  method  can  be 
consistent  and  a useful  adjunct  for 
identification  of  failure  origins. 
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1.  Introduction 

In  previous  work  [1,2]  bending  tests  were  con- 
ducted on  adhesively  bonded  specimens  of  a dental 
alloy.  The  purpose  was  to  determine: 

a)  how  much  the  bond  thickness  influenced  the 
test  results; 

b)  whether  the  failure  origins  appeared  to  be  the 
same  for  the  two  different  bond  thicknesses  em- 
ployed; 

c)  failure  origin  through  analysis  using  the  risk-of- 
rupture  (RR)  function. 


In  this  paper,  the  finite  element  method  was  used 
to  arrive  at  the  stress  distributions  used  for  the  RR 
analyses  that  employ  the  Weibull  risk-of-rupture 
function1  [3].  An  ancillary  purpose,  therefore,  was 
to  ascertain  how  sensitive  the  analyses  were  to  the 
fineness  of  the  finite-element  mesh  and,  hence, 
whether  the  method  can  be  applied  with  confi- 
dence to  the  analyses  conducted  for  a,  b,  and  c 
above. 


1 Today,  a risk-of-rupture  function,  defined  by  Weibull,  would 
be  recognized  as  a hazard  function. 
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2.  Materials  and  Methods 

The  failure  of  brittle  materials  is  typically 
catastrophic  and  in  many  instances  the  failure 
stresses  obtained  from  a set  of  test  specimens  fol- 
low a Weibull  distribution.  For  a homogeneous 
isotropic  material  subjected  to  a uniform  tensile 
stress,  a,  the  probability  of  failure,  P(a),  is  given 
by: 


P(a)  = 1 -exp-[5(o-/ai),u)'"],  (1) 

where  ob,u  is  a characteristic  strength  for  a speci- 
men of  unit  dimension;  m is  a shape  factor 
(Weibull  modulus);  and  5 is  a size  factor  (the  ratio 
of  the  failure  originating  dimension  to  a unit  di- 
mension of  the  same  kind)  and  represents  the  vol- 
ume, V,  area.  A,  or  other  dimension  in  which 
reside  the  flaws  from  which  the  failure  originates 

[3] .  Eq.  (1)  is  often  written  as: 

P(cr)  = 1 — exp  — [(cr/oo)”1],  (la) 

where  the  size  of  the  specimen,  8,  is  subsumed  into 
ab.  This  form  of  the  equation  is  commonly  used 
when  analyzing  test  data  and  the  effects  of  speci- 
men size  are  ignored.  It  is  clear  from  Eq.  (1)  that 
for  specimens  of  two  sizes,  8i  and  82,  with  the  same 
failure  origins  and  presenting  the  same  distribution 
of  failure  stresses  (nn  =m2),  there  will  be  two  dif- 
ferent values  of  cr0  for  Eq.  (la),  with  the  larger  size, 
call  it  Si,  leading  to  a value,  cr(U,  that  is  less  than 
Ob, 2- 

For  such  specimen  sets,  the  relation  between  the 
characteristic  strengths  calculated  by  Eq.  (la)  is 

[4] : 


(<7/OT'-!)"  [ («/<», .„y 

J 82 


d 5, 


dS, 


(4) 


Canceling  terms  in  ao,u  on  the  right  side  and  a on 
the  left  side  leads  to 


When  o;(5)  is  not  known  as  an  explicit  function, 
the  relationship  between  0-0,1  and  0-0,2  may,  in  princi- 
ple, be  approximated  by  [3]): 


^»,2  = o-o,i[(2^A5w)/(2^A52,)]1/m,  (5) 

j 

and  the  validity  of  the  approximation  must  be 
checked  by  computation.  Here  the  summations  are 
over  all  the  elements  considered  to  be  involved 
with  the  failure  (interface,  volume,  etc.)  and  the 
stresses  can  be  evaluated  by  the  finite-element 
method  of  analysis. 

Note:  The  stress  field  in  the  bond  region  is  typi- 
cally three  dimensional;  the  analyses  of  this  paper 
utilize  unidirectional  tensile  stresses  because  alter- 
ations to  the  principal  stresses  were  found  to  be 
minor  and  may  be  ignored.  We  also  note  that  for  a 
variety  of  reasons  (plasticity,  change  in  composi- 
tion, properties,  or  flaw  populations)  this  analysis 
method  may  not  apply  for  very  thin  bond  layers 
approaching  micrometers  or  less. 


o-0,2  = 0-o,,[51/52]1/m.  (2) 

For  a nonuniform  tensile  stress  field,  a more  gen- 
eral form  of  Eq.  (1)  is  necessary: 

P{ct)  = 1 -exp-[/5(o-/o-o)"'d5],  (3) 

where  the  region  of  integration  over  8 is  the  region 
critical  to  failure  (rupture)  and  it  can  be  in  one, 
two,  or  three  dimensions.  Then  the  relationship  be- 
tween the  values  of  0-0  [Eq.  (la)]  as  determined 
from  experiments  on  sets  of  specimens  having  ei- 
ther one  or  the  other  of  the  bond  thicknesses,  is 
given  by  the  ratios  of  the  exponents  of  Eq.  (3),  i.e., 


Each  assumed  failure  origin  for  a specimen  has 
its  own  specific  8 with  its  associated  stresses.  When 
the  ratios  of  volumes,  surface  areas,  interface  ar- 
eas, edge  lengths  etc.  (any  dimensions  containing 
the  flaws  from  which  failures  may  originate)  are 
properly  chosen  to  be  different  for  experimental 
tests,  only  one  set  of  8,’s,  0/ s,  S/s,  and  oj’s  should 
produce  coincidence  between  the  experimentally 
determined  ratio  of  Ob’s  and  the  ratios  of  either  the 
integrals  shown  by  Eq.  (4a)  or  the  summations  as 
shown  in  Eq.  (5). 

As  with  any  analysis  employing  the  finite  ele- 
ment (FE)  method  for  determination  of  the 
stresses,  a critical  question  arises  as  to  the  FE- 
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mesh  sensitivity  of  Eq.  (4).  If  sensitive,  then  the 
method  would  not,  in  actuality,  be  useful  for  the 
correlation  of  results  from  differently  sized  speci- 
mens. 

To  provide  insight  into  the  ability  of  this  ap- 
proach to  identify  sources  of  failure,  rectangular 
bond  specimens  as  shown  in  Fig.  1 were  prepared 
for  testing  in  four-point  bending,  with  either  50  p.m 
or  200  pm  bond-thicknesses.  Rectangular  bars  of 
the  bulk  bonding  material  were  also  tested  in 
three-point  bending.  The  details  of  specimen 
preparation  were  given  in  a presentation  by  Keeny 
et  al.  [1].  The  number  of  specimens  and  the  results 
for  each  test  series  are  shown  in  Table  1. 

A three-dimensional,  finite-element  elasticity 
model2  was  used  for  evaluation  of  the  stress  distri- 
bution throughout  the  volume  of  the  bond  region. 


The  original  bond  model  (Fig.  2)  consisted  of  2,197 
elements  in  one  quadrant  of  the  specimen  which 
had  three  planes  of  symmetry.  Subsequent  refine- 
ments of  interface  and  surface  elements  led  to  ele- 
ments 1/2  and  1/3  the  original  size.  The  validity  of 
the  model  was  checked  by  comparison  of  the  finite- 
element  results  for  a homogeneous  beam  with  the 
analytical  solution.  Examples  of  how  the  bending 
tensile  stresses  change  as  a function  of  the  thick- 
ness of  the  bond  layer  are  shown  in  Figs.  3,  4,  and 
5,  for  which  Eb=EJ5Q  where:  Eb  is  Young’s  mod- 
ulus for  the  bonding  material;  and  £a  is  Young’s 
modulus  for  the  alloy. 

If  the  failure  stresses  are  referenced  to  the 
stresses  along  the  surface,  the  ratio  of  the  operative 
(effective)  dimensions,  8,  from  which  the  failures 
originate  are  given  by  Eq.  (6): 


Fig.  1.  Specimen  used  for  bond  testing.  The  shaded  area  repre- 
sents the  bonding  material  between  rectangular  beams  of  alloy 
that  were  bonded  together.  Two  small  projections  were  used  to 
control  the  width  of  the  bond  at  either  50  pm  or  200  pm. 


fieff.l 

5eff,2 


( 

\CTr.1 


So  AS,, 

_i 

Ioi',A82j  ’ 

i 


(6) 


where  crR  denotes  the  reference  stress. 

For  these  calculations  the  bending  stress  at  the 
surface  was  used  as  the  reference  stress  for  calcula- 
tions of  bending  strength  and  <xR,2  = crR,i.  Then, 
from  Eq.  (2)  and  Eq.  (6) 


00,2/00,1  — (Seff,l/Seff,2)  ~ 


E<AS, 

_1 

LS0/2"AS2,,J 


(7) 


which  is  equivalent  to  Eq.  (5). 


Table  1.  Uncorrected  Weibull  parameters  and  mean  strengths 


Test 

Gap  (pm) 

/Va 

0,1  (MPa)  {range  95%}b 

m {range  95%} 

4-Pt  Bend 

50 

25 

110  {106-115} 

11.3  {8.2-14.0} 

n 

200 

25 

107  {103-111} 

11.5  {8.3-14.2} 

3-Pt  Bend 

Bars  of 
bulk  bond 
material 

54 

85  {81-89} 

6.8  {5.5-7.93} 

a N = the  number  of  specimens. 

b{}  = the  associated  confidence  bounds  on  cru  and  m as  determined  from  the  data. 


2 Developed  at  NIST. 
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Fig.  2.  A 1/8  section  of  a three-dimensional  model  for  finite  element  calcula- 
tions of  stress.  The  specimens  (Fig.  1)  had  three  planes  of  symmetry. 


Fig.  3.  Result  of  finite  element  calculation  of  the  near-surface 
bending  stress  at  the  bond  midplane  that  bisects  the  bend  speci- 
mens into  symmetrical  halves. 

3.  Results 

By  the  use  of  the  right-hand  side  of  Eq.  (7)  and 
the  finite-element-derived  stress  distributions, 
characteristic  strength  ratios  were  calculated  for 
four  potential  regions  (Table  2)  where  the  failure 
of  the  bond  could  originate,  i.e.,  volume,  surface, 
interface,  and  interface-line-junction  failure 
origins.  These  were  then  compared  with  the  results 
obtained  from  the  left  hand  side  of  Eq.  (7).  For 
these  calculations,  a value  of  m=11.5  was  used 


NORMALIZED  THICKNESS 

Fig.  4.  Surface  bending  stress  across  the  surface  of  the  bond 
material,  from  one  interface  junction  to  the  opposite  one,  as 
calculated  for  three  thicknesses.  The  5 p.m  thickness  is  pre- 
sented to  illustrate  the  trend  toward  beam  stress  calculations  as 
the  thickness  approaches  zero.  The  deviations  from  beam  the- 
ory calculations  are  appreciable  for  thick  bond  specimens,  show- 
ing the  need  to  use  the  more  robust  finite  element  method  for 
the  failure  analysis  employing  the  RR  function. 

[in  Table  1,  m was  obtained  from  Newton-Raphson 
iteration  for  fitting  experimental  data  to  Eq.  (1)]. 

The  ratio  of  the  experimental  characteristic 
strengths  has  a 95%  confidence  range  of  0.955  to 
1.11.  When  this  ratio  is  compared  with  ratios  calcu- 
lated from  the  finite  element  analyses,  the  inter- 
face (Table  2)  is  identified  as  the  origin  of  failures, 
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THICKNESS  (MICROMETERS) 


Fig.  5.  Finite  element  calculated  stress  for  the  adherend-adher- 
ent  interface  surface  junction  line,  illustrating  the  dramatic  ef- 
fect of  bond  thickness  on  interface  stresses. 


Table  2.  Ratios  of  characteristic  strengths  calculated  from  risks 
of  rupture:  (strength,  200  p.m)/(strength,  50  pm).  Four-Point- 
bending 


Assumed 

failure 

origin 

Coarse 

mesh 

Refined 
mesh  1 

Refined 
mesh  2 

Volume 

1.12 

1.12 

1.12 

Interface 

.995 

1.01 

1.06 

Surface 

1.26 

1.36 

1.32 

Interface- 
surface 
junction  line 

1.24 

1.25 

1.41 

Interface- 
surface 
junction  line 
(extrapolated)3 

1.38 

1.37 

1.37 

3 Obtained  from 

extrapolations 

of  finite 

element  derived 

stresses  at  centroids  to  the  interface  between  alloy  and  bond 
layer. 

with  all  other  failure  origins  excluded.  The  row 
with  the  next  closest  match  of  strength  ratios  is  that 
for  volume  failures  which,  with  a ratio  of  1.12,  lies 
just  outside  the  95%  confidence  range,  so  this  ar- 
gument, by  itself,  is  somewhat  unconvincing.  How- 
ever, volume  failures  are  ruled  out  because  the  m 
value  of  11.5  from  the  bond  tests  differs,  at  the 
90%  confidence  level,  from  the  value  of  m = 6.8 
which  was  obtained  from  the  bulk  specimen  test 


data.  The  bulk  specimens  can  fail  only  by  volume 
or  surface  failures.  The  strength  ratio  calculations 
for  bond  specimens  rule  out  surface  failures.  The 
m value  differences  then  are  used  to  rule  out  vol- 
ume failures. 

Hence,  the  most  reasonable  explanation  is  that 
the  bond  specimens  fail  by  interfacial  failures.  This 
is  consistent  with  features  of  the  failed  specimens 
which  always  presented  regions  showing  interfacial 
debonding. 

There  is  some  possibility  that  the  strength  of  the 
bond  itself  would  depend  on  the  bond  thickness 
due  to  a change  in  material  response  (i.e.,  forma- 
tion of  plasticity).  Such  effects  obviously  cannot  be 
dealt  with  by  the  linear  elastic  analysis  presented 
and  within  the  confines  of  this  analysis,  interfa- 
cially  initiated  failure  is  concluded. 

4.  Summary 

An  analysis  by  Weibull  RR  for  bonded  speci- 
mens of  two  different  sizes  tested  in  bending  has 
shown: 

1)  Correlations  between  characteristic  strengths, 
cr0’s,  were  possible  through  the  use  of  finite-ele- 
ment derived  stresses  in  the  RR  analysis. 

2)  The  correlations  were  not  sensitive  to  the  par- 
ticular mesh  chosen. 

3)  For  the  interface,  surface,  and  interface-surface 
junction  line,  the  stress  ratio  calculations  em- 
ploying the  element  centroid  stresses  are  not 
significantly  changed  by  use  of  stresses  from  ex- 
trapolations to  the  interface.  The  largest  differ- 
ence is  for  the  interface-surface  junction  line 
and  these  are  shown  in  Table  2. 

4)  Because  the  absolute  magnitude  of  each  RR 
calculation  changes,  mesh  of  the  same  size  and 
configuration  must  be  used  for  each  set  of  com- 
parisons. 

5)  The  origins  of  failures  can  be  determined  by 
suitable  testing  and  analysis  of  different  size 
bond  specimens  and  bulk  specimens  of  the 
bonding  material.  This  involves  the  use  of  a 
combined  approach,  analysis  of  the  Ob’s  and  m 
values. 

6)  The  determination  of  failure  origin  by  this  ap- 
proach can  be  useful  for  focusing  attention  on 
the  proper  parameters  if  improvements  in  sys- 
tem strength  or  performance  are  sought. 
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